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We have measured several branching ratios for ip' decay using the data collected by FNAh 
E835 experiment during year 2000, obtaining B(ip' -» e+e") = 0.0068 ± 0.0001 ± 0.0004, 
B(ip' -» J/ipw + w~) = 0.292 ± 0.005 ± 0.018, B(tp' -» J/ipir°w°) = 0.167 ± 0.005 ± 0.014 and 
B(ip' —> J/iprj) = 0.028 ± 0.002 ± 0.002. We also present a measurement of the dipion mass distri- 
bution in the decays ip' — > J/ipmr. 

PACS numbers: 13.75.Cs, 14.40.Gx 



I. INTRODUCTION 

Branching ratios of ip' decay have been measured by 
many experiments in e + e~ collisions. The branching ra- 
tios to e + e~ and to J/ip inclusive states have been mea- 
sured also by the Fermilab experiments E760 [| and E835 
0, which studied charmonium spectroscopy in pp anni- 
hilation. 

The strong ip' decays to J/ip proceed mainly through 
the emission of soft gluons by the cc pair, with their 
subsequent hadronization Q. The low gluon momen- 
tum makes this process nonperturbative. However some 
features of the hadronic process ip' — ► J/tpX can be pre- 
dicted using conservation laws. A prediction based on 
isospin conservation is that B(ip' — > J/ipir°Tr Q )/B(ip' — » 
J/ipir + Tr~) = i. While previous measurements have al- 
ways yielded slightly larger values, further checks of this 
ratio are needed. The multipole expansion of the gluon 
field has been used to predict a value in the range 0.10 
- 0.14 for r(V>' -> J/iprf)/T{iP' -> J^tt+tt") H and 
the 7T7T invariant mass distribution 



II. EXPERIMENTAL METHOD 

E835 studies charmonium spectroscopy in pp annihi- 
lation. All cc states can be directly formed with this 
technique. Since the cross section for charmonium for- 
mation is at best ~ 5 orders of magnitude smaller than 
the total pp cross section, in order to improve the signal 
to noise ratio the charmonium signal is extracted from 
the hadronic background by detecting electromagnetic 



final states. 

if)' decays are studied here by selecting events with a 
high invariant mass e + e~ pair in the final state, coming 
from the reactions: 

ip' — > e + e~, 

ip' — > J/ipir + ir~ — > e + e~7r + 7r~, 

ip' -> J/iPn°n a -> e + e"47, 

ip' —t J/iprj — > e + e _ 27. 

The number Na of events observed for ip' decay to a 
given final state A is given by: 



N A =£■ 



0~bkg 



e A ■ J G(E) ■ a BW (E) ■ B(iP' - A)dE S j , 



(1) 

where C is the integrated luminosity, abkg the background 
cross section, €a the overall detection efficiency for the 
channel, cjbw the Breit-Wigner cross section for pp — » ip' ', 
G(E) the center of mass energy distribution and B(ip' — > 
A) the branching ratio for the decay ip' — > A. 

Measuring the ratio of branching ratios for two chan- 
nels A and B on the same data sample, several factors in 
JTJ cancel, leading to: 



B{iP' -> A) 
B{ip' -> B) 



Na - N bkg ^A 



N f 



Ni 



bkg,B 



(2) 



where N — Nbkg is the number of events after the sub- 
traction of the background contribution. The reference 
channel has been chosen to be the J /ip inclusive decay 
because it is the one for which the ratio (J2J presents the 
lowest systematics, as will be shown in Section IVT1 



2 



III. EXPERIMENTAL APPARATUS 

E835 is a fixed target experiment in which the p 
beam circulating in the Fermilab Antiproton Accumu- 
lator (AA) crosses an internal hydrogen gas jet target. 
During a typical data taking period the p are accumu- 
lated until a current of ~ 70 mA is reached. Then ac- 
cumulation stops and the p are stochastically cooled and 
decelerated to the energy of the resonance to be studied. 
At this point the H2 target is turned on and the data tak- 
ing proceeds on the resultant "stack" until the p beam 
current is around 10 mA. In this current range it is pos- 
sible to take data at constant luminosity by regulating 
the jet target density. 

After deceleration, the mean center of mass energy E cm 
of the pp system is known with a precision better than 
~ 100 keV, where the E cm distribution is gaussian with a 
a typically ranging from 200 to 400 keV. Further details 
about the AA operation during E835 data taking can be 
found in [To| . 

The E835 detector is a nonmagnetic spectrometer with 
cylindrical symmetry around the beam axis. The inner 
part of the detector is the charged tracking system; for 
the year 2000 data taking it was composed of three cylin- 
drical hodoscopes, two straw chambers for the measure- 
ment of the azimuthal angle cp (around the beam axis), 
two scintillating fiber detectors for the measurement of 
the polar angle 9 (with respect to the beam axis) and 
an additional hodoscope in the forward direction used 
as a charged veto. The three cylindrical hodoscopes are 
segmented in <p in 8 (for the inner), 24 (for the interme- 
diate) and 32 (for the outer) modules respectively. The 
hodoscopes provide dE/dx information for the charged 
tracks and they are used to form the charged hardware 
trigger. 

A 16 cell Cerenkov counter with 8 azimuthal and two 
polar angle sections cover the full azimuth and the polar 
angle region 15° < 9 < 65° and allows the separation of 
high energy e , mostly produced in J/ip and ip' decay, 
from the other charged particles, which is used in the 
first level trigger. 

Two electromagnetic calorimeters, the Central 
Calorimeter (CCAL) and the Forward Calorime- 
ter (FCAL) cover the region 11° < 9 < 70° and 
3° < 9 < 11° respectively. The CCAL energy resolution 
is <je/E = 1.4% + 6%/y/E{GeV) and its angular 
resolutions are er^ ~ 11 mrad and ag ~ 6 mrad. FCAL 
is not used in the following analysis. 

The luminosity is measured with ~ 2% accuracy by 
means of solid state detectors located below the interac- 
tion region at 9 ~ 90° to the p beam. 

IV. EVENT SELECTION 

In year 2000, E835 collected a total integrated lumi- 
nosity of L = 113 pb _1 of data. The data sample used 
in this analysis consists of 14.4 pb _1 collected in the 



ip' resonance region, with 12.5 pb" 1 on the resonance and 
1.9 pb _1 in two samples at energies above (3704.9 McV) 
and below (3666.1 MeV) the ip' , used for background 
measurements. 

Interesting events are characterized by two nearly 
back-to-back high energy e^. The first level hardware 
trigger is the logical OR of one main trigger condition 
and two control triggers (ll). The main hardware trig- 
ger requires two "electron tracks" as defined by signals 
in the Cerenkov and corresponding hodoscope elements 
in coincidence with two back to back energy deposits in 
the calorimeter. The efficiency of the main trigger was 
measured on data taken with a dedicated trigger, and 
continuously monitored by the two control triggers, the 
first with no Cerenkov requirement and the second with 
no calorimeter requirement. 

As a preliminary selection, all events with 
e + e~ candidates with invariant mass m ee < 2.6 
GeV are rejected. 

A maximum likelihood method called "Electron 
Weight" (EW) has been developed for the rejection 
of backgrounds that mimic e tracks in the detector 
(mainly 7 conversion and 7r° Dalitz decay), comparing 
the signal and background probability for the candidate 
e + e~ pair on the basis of pulse heights in hodoscopes and 
Cerenkov counters and the CCAL shower shape. Figure 
n shows the logarithm of the EW product for the two e 
candidates for data in the ip' energy region and for the 
background. We choose to cut at EW\ ■ EW2 > 1.5. In 
this way a clean sample of e+e~ due to charmonium de- 
cay is selected. Figure |21 shows the invariant mass distri- 
butions of the e + e _ candidates before and after the EW 
cut. After the cut, the two m ee peaks due to the direct 
decay ip' — ► e + e _ and the cascade ip' — + J/ipX — > e + e~ X 
are clearly visible. 

All the events containing an additional shower in the 
CCAL within 100 mrad of the e ± tracks are rejected. 
This cut yields homogeneous values of the efficiency cew 
for all the final states and reduces the systematic error 
in the final branching ratio measurements. 

ip' — * e + e _ and ip' — > J/ipX are first selected using 
kinematical fits requiring prob(x 2 ) > 10~ 4 ; the ip' — > 
J/ipX are then tagged as J/ipir + n~, J/i/jir°7r° or J/iprj 
with kinematical fit selections using the information on 
the additional tracks in the detector. 

For J/ipn + 'K~ the ^ tracks are reconstructed within 
the charged tracking system acceptance (15° < 6„± < 
55°) combining the (p measurements obtained with the 
straw chambers (0-lines) and the 9 measurements ob- 
tained with the scintillating fiber detectors (0-lines). For 
each event containing at least two <p- and two 0-lines in 
addition to e + e~ tracks, all the possible pairs of charged 
9-<p lines combinations are tested with the kinematical 
fit as 7r + 7r~ candidates. The combination with the high- 
est x 2 probability is chosen. The ^ energy cannot be 
measured in the E835 detector so a 3C kinematical fit is 
applied for the selection of J /ipir + ir~ events. 

The 71"^ have relatively low momentum in the labo- 
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FIG. 1: log 1Q (EWi ■ EW2) distribution for events with m ee > 
2.6 GeV in the ip' energy region (open) and on the back- 
ground with normalized luminosity (cross hatched). The cut 
log 10 (£Wi • EW2) > log 10 (1.5) is indicated by the arrow. 



FIG. 2: m ee distribution for all the events with candidate 
e + e~ pair with m ee > 2.6 GeV (open). The events se- 
lected with the EW\ ■ EW2 > 1.5 cut in the tp' energy region 
(hatched) and on the background scaled with the luminosity 
(cross hatched), are also shown. 



ratory and the multiple scattering effect dominates the 
uncertainty in the measurement of their directions. The 
associated error has been determined using a full detector 
simulation based on the GEANT 0] package. 

Photons are detected as energy deposits in the CCAL 
fiducial region 12° < 9 7 < 68°. ir° and r\ candidates 
are sought in their 77 decay, looking for photon pairs 
of invariant mass between 50 MeV and 350 MeV or 300 
MeV and 800 MeV respectively. 

The previous final state categories are mutually exclu- 
sive. If an event passes more than one of the previous 
selections, the event is assigned to the final state with 
the highest prob(x 2 ). The number of events selected is 
reported in Table [I] We find that of the 31193 candi- 
date events, 30492 (98%) are selected as e + e~ or J/tpX, 
in good agreement with the Monte Carlo results for the 
efficiencies of the two selections, as described in Section 

EH 



pared to the resonance signal in Figure |21 The number 
of background events from this source to be subtracted 
is obtained scaling the number of events observed in each 
sample with the integrated luminosity. 

The misidentification background is due to resonant 
events which are not correctly classified. It is evaluated 
using a Monte Carlo simulation: a sample of events with 
full detector simulation is generated for all the contam- 
inating channels and analyzed with the same analysis 
used for the data. The probability of event misidenti- 
fication is shown in Table [H] The number of internal 
background events to be subtracted for each exclusive 
contaminating channel (Ni nt (A), where A indicates the 
contaminating final state) is obtained from the number 
of observed J/i/jX events, scaled with the branching ratio 
decay mode obtained from the PDG 0] , and multiplied 
by the misidentification probability P m isid (Table ITT|) : 



V. BACKGROUND SUBTRACTION 

Two sources of background must be subtracted from 
the observed numbers of events: a non resonant back- 
ground contribution and an internal background due to 
final state misidentification. 

The non resonant (external) background is due to 
events not coming from ip' decay. It is measured on data 
taken in the energy regions far from the resonance (stacks 
54a and 54c shown in Table P) by applying the same 
analysis as for the ip 1 peak data; its contribution is com- 



N int (A) = N(J/^X) 



B{f -» A) 



B(i/j' -» J/ipX) 



PDG 



(3) 

The numbers of events after background subtraction 
are summarized in Table UTT1 



VI. SELECTION EFFICIENCY 

The overall detection efficiency e for each exclusive 
channel is the product of the efficiencies of the cuts used 
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TABLE I: Number of events observed in each decay channel in E835 year 2000 if)' data, divided by stack. Stacks 54a and 54c 
are collected far from the resonance and are used for the background measurements. The charged tracking system was off for 
part of the stack 54a. During this period the J/tpir + ir~ could not be detected so the effective total integrated luminosity on 
the background (stacks 54a and 54c) for this channel is 1.7 pb _1 . 



Stack 


C (pb" 1 ) 


Candidates 


e + e 


J/ipX 




J/^7r°7r° 


J/lpTj 


1 


0.748 


831 


155 


658 


114 


31 


12 


2 


1.008 


3595 


636 


2893 


511 


187 


41 


14 


0.992 


2248 


388 


1807 


323 


132 


27 


29 


0.992 


2083 


313 


1706 


273 


126 


17 


30 


0.396 


951 


158 


769 


135 


45 


9 


49 


2.566 


5890 


931 


4840 


781 


322 


67 


50 


1.275 


3633 


597 


2966 


385 


205 


43 


51 


2.103 


5016 


812 


4080 


660 


275 


65 


54b 


2.401 


6946 


1114 


5669 


971 


373 


82 


54a 


1.153 


35 


1 


29 


1 


1 





54c 


0.780 


23 





19 


1 


1 






TABLE II: Misidentification probabilities. Contaminations smaller than 10 3 have been neglected. 



Generated 


e+e- 




J/tpX 


Tagged as 
J/ipir + ir~ 




J/iprj 


e + e~ 






0.038 ± 0.001 


< 10" 3 


< 10" 3 


< 10" 3 


J /xjjn + iT~ — ► e + e~ ix + ix~ 


< 10" 


3 






< 10" 3 


< 10~ 3 


J/ipTV°Tr° — > e + e~7777 


< 10" 


3 




0.0015 ± 0.0003 




0.0056 ± 0.0005 


J/iprj — > e + e~77 


< 10" 


3 




< 10~ 3 


< 10~ 3 




J/ipr] — » e + e - 7T + 7r - 7r° 


< 10" 


3 




0.104 ±0.002 


< 10 -3 


< 10" 3 


J/tpT] -» e+e-TrW 


< 10" 


3 




0.0031 ± 0.0004 


0.061 ± 0.002 


< 10~ 3 


Xci7 — * J/i>11 — * e + e"77 


< 10" 


3 




< 10" 3 


< 10" 3 


0.0150 ± 0.0008 


Xc2j — > J/ iA77 — > e + e~77 


< 10" 


3 




< 10~ 3 


< 10~ 3 


< 10~ 3 



for the event selection: 

e = a ■ e tr ig ■ e mee ■ cew ■ £sel, (4) 

where a is the acceptance for the e + e~ pair coming from 
charmonium decay, etrig the trigger efficiency, e mee the ef- 
ficiency of the e + e~ invariant mass cut (M ee > 2.6 GcV), 
e E w the efficiency of the EW cut (EW X ■ EW 2 > 1.5) 
and e se i the final selection efficiency, which includes also 
the acceptance for all the remaining particles of the final 
state. The relevant quantity in this analysis is the ratio 
of the efficiencies for different channels. 



A. Acceptance and trigger efficiency 

The final state is characterized by a high invariant mass 
e + e~ pair. The hardware e + e~ trigger is based on the 
Cerenkov counter, multiplicity and topology in the ho- 
doscopes and multiplicity and topology and energy re- 
lease in the CCAL. 

The e + e~ acceptance is determined by the Cerenkov 
fiducial volume, which is defined as 15° < 9 < 60°. 

In this analysis, the product a ■ etrig is evaluated for 
each channel using a sample of Monte Carlo events. 

The main source of error on a ■ e tr i g is the uncertainty 
in the e + e _ angular distribution. In the reaction pp — > 



ijj' — > e + e~ the e + e~ pair is distributed according to 1 + 
Xj/)i cos 9* , where 9* is the angle between the e^- and the 
beam directions in the center of mass system and A^,' 
is the angular distribution parameter. In what follows 
we used = 0.67 ± 0.16, as recently measured by this 
experiment |14| . 

For the decay ip' — > J/ipirir, the transition is dominated 
by the emission of an L — dipion . In this case the 
J I tp has the same polarization as the tp' and the angular 
distribution for the J/ i/j decay is 1 + cos 9* ; in the 
Monte Carlo simulated events we also assumed S wave 
between the dipion and the J/tjj . 

For the J/iprj channel, the e + e~ angular distribution 
is 1 + |A^' — X^f cos 91 as discussed in [T^ |. 

For the double radiative decay channels ip' — > Xdl — > 
J/^77, representing a background for some of the chan- 
nels analyzed here, the ip' radiative decay angular distri- 
butions are simulated using the quadrupole amplitudes 
measured by CBAL [13, while for the Xcj radiative de- 
cays, the angular distributions used in the simulation are 
the ones measured by E835 [Isj . 

The results obtained for a ■ e tr i g are summarized in 
Table IIVI Using these values we obtained the ratios of 
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TABLE III: Numbers of events with statistical errors and background subtraction. N ex t and Ni n t are the external and internal 
background contributions. 

e+e~ J/iPX J/Tpiv+TT- J/^ttV' J/t/jt? 

N evta 5104 ± 71 25388 ± 159 4153 ± 64 1696 ± 41 363 ± 19 

N ext 7tf 316 ±46 15±ig 13^ 7 < 11 (68% C.L.) 

N int - 206 ±9 48 ±4 28 ± 2 65 ± 7 

N 5097 ±73 24866 ±166 4090 ± 67 1655 ± 44 298 ± 20 



TABLE IV: a ■ e tr i g for different values of X^i obtained with 
Monte Carlo simulation. The Xc./7 channel, by radiative de- 
cay, leads to 7/7/177. The values for the J ftp inclusive decay 
are obtained as the weighted average of the exclusive chan- 
nels. 



Channel 


X r = 0.51 


V = 0.67 


X^i = 0.84 


e+e" 


0.5114 


0.4995 


0.4893 


J/ipn + ir~ 


0.4820 


0.4681 


0.4577 


J/i>TT°7T° 


0.4831 


0.4738 


0.4619 


J/ipr) 


0.5672 


0.5751 


0.5823 


XcOl 


0.5249 


0.5248 


0.5233 


Xcll 


0.5094 


0.5068 


0.5030 


Xc27 


0.5026 


0.4937 


0.4923 


J/ipX 


0.4890 


0.4785 


0.4692 



the a ■ ttrig product between different channels, namely: 



a ■ e trig (J/i>X) 
a ■ e tri g(e + e-) 
a ■ etrigjJ/ipX) 
a ■ e trig (J/ipTT+Tr-) 
a ■ e t „ g (J/ipX) 

a ■ etrtgiJ/lpTT ^ ) 
a ■ etrigjJ/lpX) 
& ' ^trig 



0.9579 ±0.0016, (5) 

1.022 ± 0.007, (6) 

1.010 ±0.006, (7) 

0.83 ±0.03, (8) 



where the error is systematic and comes from the uncer- 
tainty in the e + e~ angular distributions. 



B. e e invariant mass cut efficiency 

The efficiency of the e + e~ invariant mass cut, m ee > 
2.6 GeV, is e mee ~ 96%, and affects all the channels 
in the same way, except ip' — > e + e~ which has a higher 
efficiency because of the higher e + e~ invariant mass. The 
efficiency ratio is determined from Monte Carlo to be 



.(j}'->J/ipX->e + e~X) 
e (f- 



0.992 ±0.001. 



C. Electron weight cut efficiency 

The presence of additional showers in the CCAL in 
proximity to the e ± tracks modifies the CCAL shower 
shape for the e^, resulting in a lower EW cut efficiency 



for these events. The events containing extra CCAL 
showers within 100 mrad of tracks are rejected; this 
allows the use of the same EW efficiency for all channels. 
The residual fluctuations of the cew values for differ- 
ent channels, due to additional effects, have been studied 
with the Monte Carlo and are within 2%, which has been 
taken as the EW cut systematic error. 



D. Final selection efficiency 

The effect of the 100 mrad cut in the CCAL, the accep- 
tance for exclusive reactions, and the kinematic fit selec- 
tion efficiency for each analyzed channel are included in 
the final selection efficiency e se i . This is determined from 
Monte Carlo event samples by performing the same cuts 
applied on data. The systematic error associated with 
the efficiency is obtained by trying various prob(% 2 ) cuts 
in the range 10~ 4 - 10 _1 for the J/iprj channel and 10 -6 - 
10~ 2 for all the other channels, for both data and Monte 
Carlo events, and observing the result fluctuations. The 
results obtained for each channel are shown in Table [V] 



VII. RESULTS 

Using the efficiency values obtained in the previous 
section and (0) we obtain: 



B(ip' -> WO ■ B(J/ip e+e-) 
= 0.206 ±0.003 ±0.008, 

B(ip' J/ipX) 
= 0.525 ±0.009 ± 0.022, 
B{ji' -> J/^°7r°) • B(tt° -> 77 ) 2 
B(ip' -» J/tpX) 
= 0.292 ±0.008 ±0.019, 
B(ip' -> J/iprj) ■ B(r) -> 77) 

B(ip' ^ J/ipX) 
= 0.0197 ±0.0013 ±0.0013, 



(9) 
(10) 

(11) 
(12) 



where the first error is statistical and the second system- 
atic. 

The number of events observed in the J/ipir°ir° and 
J I iprj decay modes could also be normalized to the events 
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TABLE V: e 3e i for each decay channel. The error is systematic and is obtained by examining several prob(x 2 ) cuts. 



Channel 


e+e 


J/ipX 


J/ipir + ir 




J/tprj 


(■se.1 


0.943 ± 0.031 


0.996 ±0.011 


0.319 ±0.011 


0.229 ±0.014 


0.502 ± 0.024 



observed in the J/ifm + ir . Our result for these decay 
modes could alternatively be expressed as: 



B(tp' -> J/iPtt°tt°) 



0.571 ±0.018 ±0.044, (13) 



to be compared with the value 0.5 expected from isospin 
conservation, and: 



0.095 ±0.007 ±0.007, (14) 



B(tp' -*■ J/jjrj) 
B^f? ^ J/tJjtt+t:-) 



in agreement with theoretical estimates 

From our measurement and the world averages for 
B(V>' -> J/ipX) = 0.557 ± 0.026, B{J/ip -> e+e") = 
0.0593 ± 0.0010, B(tt° -> 77) = 0.98798 ± 0.00032 and 
B(rj -> 77) = 0.3943 ± 0.0026 we can derive the 
following values for the ip' branching ratios: 



B(if>' e+e-) 
B(V/ -> J/^7r°7T ) 



(6.8 ±0.1 ±0.4) x 10~ 3 ,(15) 
0.292 ±0.005 ±0.018, (16) 
0.167 ±0.005 ±0.014, (17) 
0.028 ±0.002 ± 0.002. (18) 



The results are in excellent agreement with the recent 
measurements published by BES E3- The present 
result is compared with the ones obtained previously by 
our experiment and E760 in Table IVII The larger sys- 
tematic error in the present results comes from a more 
conservative evaluation of the kinematic fit systematic. 



VIII. DIPION INVARIANT MASS 
DISTRIBUTION MEASUREMENT 

The dipion invariant mass {m^) distributions for 
J/ipiT + TT~ and J/ , 07r O 7r events, corrected for the detec- 
tor acceptance, are shown in Figure 

A possible parametrization for the m„ distribution is 
0: 



dT 

dm^. 



2 \2 



oc Ps ■ (m^ - A^mJ 



(19) 



where the phase space P$ is: 



2 

7T7T 



The measured 771^+^- and m^o^o distributions, cor- 
rected for the detector acceptance, are fitted to the func- 
tion i|19|) convoluted with their resolution as determined 
by Monte Carlo. 

The fitting procedure has been verified to reproduce 
the correct (input) result when applied to simulated 
events, and to perform equally well on the neutral and 
charged dipion modes. 

The results obtained for the distribution are 

A OT+7r - = 3.31 ± 0.15±g;?l {x 2 /ndf = 32.3/29 = 1.1) 
and A^o = 4.06 ± 0.25±S;?| (x 2 /ndf = MA/29 = 1.9) 
where the first error is statistical and the second system- 
atic. The m„ distributions for J/ipir + ir~ and J/ipTr°ir° 
are expected to be the same. 

For the evaluation of the systematic error, several fits 
to the data have been done, varying the binning, the 
parameters of the resolution function and the acceptance 
correction factors. The value of A^ was found to be very 



sensitive to the parameters of the resolution function, 
which is described by a double gaussian whose parame- 
ters are obtained from the Monte Carlo. The systematic 
errors have been obtained by varying the parameters of 
the distribution to take into account possible discrepan- 
cies between the data and the Monte Carlo. For the 
acceptance correction factors, some sets of values have 
been calculated under the hypothesis of A^ values from 
2 to 6. Then the data have been fitted using different cor- 
rections; the Ajr,,- values fitted on data using the different 
correction factors show variations by a fraction of the 
statistical error, and so no systematic error is associated 
with the acceptance correction. We also observe that the 
systematic uncertainty introduced by A Wff on the accep- 
tance used in the measurement of the branching ratios of 
the J I tpTTTi channels is negligible with respect to the ones 
in Tables llVl and M 

The only other available measurement of A nn has been 
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TABLE VI: V' BR measurements obtained by E835. 





B(ip' -> e+e-) 


B(tp' -> J/iPir+n-) 


->• J/ip-K°-R°) 


B{i>' -» 




B(ip' -> J/tpX) 


B(ip' -> J/ipX) 


B(ip' -» J/tpX) 


B(ip' -> J/ipX) 


E760 [1] 


0.0144 ± 0.0008 


0.496 ± 0.037 


0.323 ± 0.033 


0.061 ± 0.009 


E835 [2j 


0.0128 ± 0.0004 




0.328 ±0.015 


0.072 ± 0.009 


E835 (this paper) 


0.0122 ± 0.0002 ± 0.0005 


0.525 ± 0.009 ± 0.022 


0.300 ± 0.008 ± 0.022 


0.050 ± 0.006 ± 0.003 




FIG. 3: 7T7T invariant mass distribution for (a) J/ipn+n' events, corrected for the acceptance. The solid line is 

the fit with Ps ■ (m 2 ^ — A^m^) 2 . 



obtained by BES for the decay tp' — » J/ip-K+ir , which 
yielded A 7r+7r - = 4.35 ± 0.06 ± 0.17 15]. The fit to our 
data using this value yields a \ 2 /ndf = 2.7. 
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